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The possible existence of stable axial octupole and tetrahedral deformations is investi- 
gated in 80 Zr and 98 Zr. HFBCS calculations with parity projection have been performed 
for various parametrizations of the Skyrme energy functional. The correlation and ex- 
citation energies of negative parity states associated with shape fluctuations have been 
obtained using the generator coordinate method (GCM). The results indicate that in 
these nuclei both the axial octupole and tetrahedral deformations are of dynamic char- 
acter and possess similar characteristics. Various Skyrme forces give consistent results 
as a function of these two octupole degrees of freedom both at the mean-field level as 
well as for configuration mixing calculations. 



1. Introduction 

It has been recently conjectured that many nuclei throughout the periodic table 
possess a tetrahedral deformation in their ground- or low-lying isomeric state. This 
type of deformation is realized mainly through the non-zero intrinsic octupole mo- 
ment Q32 oc ?" 3 (Y32 + Y3-2) accompanied by vanishing quadrupole deformation. 
The first study pointing at the importance of the tetrahedral degree of freedom 
in many-fermion systems has been reported in Ref. . Recently more realistic ap- 
proaches based on the microscopic-macroscopic model with Woods-Saxon Hamilto- 
nian revealed that several nuclei may possess stable tetrahedral deformation. These 
nuclei are believed to form islands on the nuclear chart around " tetrahedral magic" 
numbers of neutrons and protons: 16,20,32,40,56 - 58,70,90 - 94 2 < 3 < 4 . The en- 
hanced susceptibility towards the tetrahedral deformation is associated with the 
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symmetry of pyramid- like shapes. It leads to the appearance of two- and four-fold 
degeneracies in the single-particle spectrum and consequently generates large shell 
effects stabilizing this deformed configuration 3 < 2 ' 4 * 5 . 

On the other hand the energy difference between the spherical and tetrahedral 
configurations is a sensitive function of the pairing strength. Consequently a stable 
tetrahedral deformation is a result of a delicate balance between shell effects and 
pairing correlations. Indeed it was shown that for two tetrahedral magic nuclei: 80 Zr 
and 98 Zr, the energy difference between spherical and tetrahedral configurations 
does not exceed 1 MeV. Consequently quantum fluctuations beyond the mean-field 
play a significant role 6 . Moreover these nuclei also exhibit other types of octupole 
deformations: the coupling between the neutron <i 5 / 2 and h\\ji orbitals and the 
proton p 3 / 2 and t/g/2 orbitals leads to both axial and non-axial octupole correlations 
7,6 . Moreover octupole deformations are in competition with the quadrupole mode. 

In the present article we analyze the influence of various parametrizations of the 
Skyrme interaction and pairing strengths on the existence of tetrahedral deforma- 
tion in 80 Zr and 98 Zr. In the first section the competition between axial octupole 
and tetrahedral shapes at the mean-field level is discussed. In the next section 
the nuclear dynamics beyond the mean-field is investigated. This study is directed 
towards the determination of the correlation energies and excitation energies of 
negative parity states. 

2. Hartree-Fock + BCS approach and the parity projection. 

The Hartree-Fock + BCS (HFBCS) method has been applied to obtain the energy 
of a nucleus as function of octupole degrees of freedom. The Hartree-Fock equations 
have been solved on a 3-dimensional mesh in coordinate space. The details of the 
calculations can be found in Refs. 8 ' 9 ' 6 . The pairing interaction has been treated in 
the BCS approximation including the Lipkin-Nogami (LN) correction 10 . A zero- 
range density-dependent pairing interaction has been used: 

V V air = \gi{\ - P«)S(r - r') (l - , (I) 

where i — n,p for neutrons and protons, respectively. As in previous applications, 
we set po = O.I6fm~ 3 . The strength of the pairing force has been adjusted to 
"experimental" pairing gaps, extracted from the odd-even mass staggering using a 
three-point filter from Ref. 11 . 

The behavior of the HFBCS energy as function of the axial octupole and 
tetrahedral degrees of freedom qualitatively agrees with the results presented in 
Ref. 12 > 13 > 14 . The existence of axial and non-axial octupole minima is a sensitive 
function of the pairing strength. The minima are not well pronounced and vanish 
when pairing is increased 6 . However the HFBCS approach does not conserve nei- 
ther parity nor particle number. Namely, the HFBCS solutions represents a mixture 
of positive and negative parity states. 
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Fig. 1. The parity-projected mean-field energy as function of shape parameters /33 M (see text 
for definition) is shown for three parametrizations of the Skyrme force. The two figures at the 
left-hand side show the energy as function of the axial octupole shape parameter fizo, whereas 
the figures at the right-hand side show the energy as function of (3^2 ■ The solid lines with filled 
symbols denote positive parity solutions and the dotted lines with open symbols denote negative 
parity solutions. Circles, triangles and squares correspond to the Skyrme parametrization: Sly4, 
SkM* and SHI, respectively. The energies are shown relative to the lowest energy of the spherical 
configuration. 



The restoration of broken symmetries has been performed through the projec- 
tion on a definite parity and particle number. The corrected energies for N neutrons, 
Z protons and both parities are obtained by exact projection: 

m^w^MwawMM, (2) 

where |<X/?3m)) are HFBCS wave functions generated with the constraint 
(<M/33m)I ( ?3mI ( /K/53m)) = Cpfop. The shape parameters /3 3/J are related to the oc- 
tupole moments through the relation: /3 30 = {Qzo)/Cq, (3^2 = (<332)/C2, where 
3 

Co = —A 2 rl, C'2 = C0/V2 with r = 1.2/m. The operator P(±,n.z) is the product 
of operators projecting on it = ±1 parity and on N neutrons and Z protons. 

The parity-projected energies for three parametrizations of the Skyrme force 
are shown in the Fig. 1. One may notice that all forces give qualitatively the same 
dependence (apart from a trivial energy shift) as function of the axial octupole 
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Fig. 2. The parity projected mean-field energy as function of shape parameters fizp, obtained using 
the Sly4 force. Two values of the pairing strength have been used: reproducing the experimental 
gaps (circles) and with pairing gaps increased by a factor 2 (squares). The solid (dotted) lines 
with filled (open) symbols denote positive (negative) parity solutions. 




and tetrahedral degree of freedom. As usual after parity restoration 15 the energy 
minima for positive parity states are shifted towards smaller octupole deformations 
compared to the HFBCS minima, while the negative parity states have larger de- 
formations. For both nuclei the energy minima for positive parity correspond to 
very similar /3 30 and ^32 values. For the negative parity curve, /3 32 is systematically 
larger than /3 30 in the minimum. 

It is instructive to investigate the sensitivity of the parity projected solutions 
to the pairing strength. In the Fig. 2 we have shown the parity projected energies 
for the SLy4 parametrization of the Skyrme force for normal and increased (twice) 
strength of pairing correlations. Note that contrary to the behavior of unprojected 
HFBCS energies, in this case the dependence of the energy as function of octupole 
degrees of freedom practically remains unchanged. In particular, the position of 
the energy minima for both negative and positive parity states remains unaltered. 
Similar results have been obtained for other Skyrme parametrizations. 
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Fig. 3. The square of the modulus of the collective wave function vs. shape parameters /33 M are 
plotted for three parametrizations of the Skyrme force. The two figures at the left-hand side show 
the energy as a function of the axial octupolc shape parameter /?3o, whereas the two at the right- 
hand side show the energy as function of /?32 (see text). The symbols used for each curve are the 
same as in Fig. 1. 

3. Generator coordinate method. 

Mean-field results suggest that shape fluctuations are important in this case. In or- 
der to quantify this effect we have applied the generator coordinate method (GCM). 
It allows to calculate the correlation energies associated with shape fluctuations and 
to determine the structure of the collective wave functions in terms of the contribut- 
ing mean-field configurations. Namely, a collective wave function is constructed by 
mixing the mean-field states corresponding to different values of the octupole mo- 
ment, after their projection on particle number and parity: 

l*> = j ' f(^)P { ±,N,z)\^))d(3 3 ^. (3) 

The coefficients /(/?3 M ) are determined by minimization of the total energy of the 
collective wave function |\I>). In practice, the integral is replaced by a discrete sum- 
mation over /3 3/J , with a number of points large enough to obtain results independent 
of the discretization 16 ' 6 . The discretized Hill- Wheeler (HW) equation was solved 
separately for each collective coordinate Q30 and Q32. 
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Fig. 4. The square of the modulus of the collective wave function vs. shape parameters are plotted 
for the Sly4 force. Two values of pairing strength have been considered: reproducing the experi- 
mental pairing gaps (circles) and with pairing gaps increased by a factor 2 (squares). The solid 
(dotted) lines with filled (open) symbols denote positive (negative) parity solutions. 

In the Fig. 3 the collective wave functions (related to f(03/j,) by an integral 
transformation) are plotted for three parametrizations of the Skyrme force. One 
can see that all Skyrme forces give consistent results. The wave functions are spread 
around the minima of the projected mean-field energy curves, with a shape typical 
of a vibration in a 1-dimcnsional energy well. Note that there is no indication that 
any particular tetrahedral configuration contribute the most to the wave function. 
The large spread of the wave function confirms the importance of shape fluctuations. 
It is reflected also in the large value of the correlation energies which are listed in 
the table 1. Correlation energies have been calculated from the prescription: 

Ecorr = E{N,Z, spher. ) - E+ , (4) 

where E(N, Z, spher.) is the energy of the particle number projected spherical con- 
figuration obtained in the HFBCS approach, and E + is the lowest positive-parity 
energy obtained in the GCM. 

Note also that the collective wave functions for both axial octupole and tetra- 
hedral coordinates have a very similar shape. Indeed the calculated correlation en- 
ergies and excitation energies of the negative parity state have very similar values, 
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Tabic 1. Results of the GCM calculations for 
three Skyrme forces: SLy4, SIII and SkM*. 
E CO rr(Q A/j) represents the correlation energy as- 
sociated with shape fluctuation described by 
multipolc moment (see eq. 4). Similarly 
EexciQxfi) denotes the excitation energy of the 
first negative-parity collective state with respect 
to the first positive-parity state. All values are in 
MeV. 



Nucleus 


SLy4 




SIII 


SkM* 


80 Zr 


E C orr(Qz()) = 


1.507 


1.459 


1.290 




E C orr(Qs2) = 


1.567 


1.499 


1.338 




E e xc(Qm) = 


2.802 


2.520 


3.111 




E exc (Q32) = 


2.425 


2.220 


2.881 




E C orr(Qw) = 


1.389 


1.387 


1.554 




E C orr(Qz2) = 


1.400 


1.485 


1.564 




Eexc(Q30) = 


2.644 


1.090 


2.116 




Eexc(Q32) = 


2.498 


0.784 


1.776 



slightly favoring the tetrahedral configuration. 

In order to check the sensitivity of the GCM results on the magnitude of pairing 
correlations, we have performed calculations with pairing strengths that produce 
gaps larger or smaller by a factor 2. Qualitatively, the GCM results are not affected 
which can be seen in the Fig. 4. The correlation energies vary by 10 — 20%. A 
doubling of the pairing gap results in approximately twice larger excitation energy 
for the negative parity states. 

4. Conclusions. 

• The existence of a stable tetrahedral deformation at the mean-field level is 
rather unlikely for 80 Zr and 98 Zr. 

• The parity projection induce a small tetrahedral (and also axial octupolc) 
deformation which is relatively independent of the pairing strength. 

• Shape fluctuations play an important role and significantly contribute to 
the correlation energy. 

• Both axial octupole and tetrahedral states have very similar characteristics 
although the tetrahedral configuration is slightly more favored. 

Summarizing we conclude that the existence of the tetrahedral deformation is of 
the dynamic character. 
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